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By using the series expansion techniques, we study the ex- 
citation spectrum for the two-dimensional quantum spin sys- 
tems with ladder, plaquette and mixed-spin structures. We 
calculate the spin excitation gap and thus determine the phase 
boundary between the spin-gap phase and the magnetically 
ordered phase rather precisely. It is found that the phase 
diagram obtained improves fairly well the one previously ob- 
tained via the ground-state susceptibility. 



I. INTRODUCTION 

Two-dimensional (2D) antiferromagnetic quantum 
spin systems with spin gap provide a new interesting 
paradigm for quantum phase transitions. A typical £x^ 
ample is the spin plaquette system such as CaV409,crEI 
which may be described by the 2D Heisenberg model with 
a plaquette structure. Another interesting example found 
recently is the 2D spin system composed of orthogonal 
spin dimers such as SrCu2(B03)2 Jijil3 which may be de- 
scribed by the 2D Heisenberg model on a square lattice 
with some diagonal exchange couplings. In these 2D spin 
systems, the plaquette or dimer structure is essential to 
stabilize the non-magnetic phase with the spin gap. Con- 
cerning the spin gap formation, the topological nature of 
spins is also important in low-dimensional systems. In 
this context, mixed-spin systems have attracted consid- 
erable attention recently, in which the spatial arrange- 
ment of different spins plays a crucial role to generate 
the spin gap or induce an antiferromagnetic long-range 
order. For the references for experiments^ 

and theoriealaila in ID cas^. 

In the previous paperjlj we have investigated the 
ground state quantities for the 2D spin systems with lad- 
der, plaqiiette-and mixed-spin structures, by extendiue 
the workgjijQtj based on the series expansion methods.EJ 
Although the quantum phase transitions have been de- 
scribed qualitatively well, it has turned out that the ob- 
tained results lead to unsatisfactory estimates for—the 
phase boundary in some region of the phase diagram.Ej A 
more crucial problem raised is to what extent our series 
expansion correctly captures the lattice structure and/or 
the spin topological properties, since our expansion ap- 
proach has relied on the lower-order expansions in cou- 
pling constants. Not only to resolve this problem but 
also to confirm our approach to be reliable, it is desirable 
to produce more accurate results by improving the series 
expansions, and also to study other quantities besides the 
ground state quantities. 



The purpose of this work is to study the excitation 
spectrum for the 2D quantum spin systems with the 
above-mentioned structures, and clarify the role of the 
competing interactions in the disordered phase with the 
spin gap. We shall see that the excitation spectrum cal- 
culated in higher orders than the ground-state suscepti- 
bility improves the phase diagram, and at the same time 
confirms that our series-expansion approach indeed pro- 
vides reliable results. 

The paper is organized as follows. In §2, we 
briefly suoimarize how to apply the series expansion 
techniquesE3 to our systems. By—performing the series 
expansion for the excited statescJ and employing the 



asymptotic analysis of power-series expansionaij in §3, 
we obtain the dispersion relations and the phase dia- 
gram. We then discuss how the spin gap phase competes 
with the magnetically ordered phase for three kinds of 2D 
quantum spin systems mentioned above. The last section 
is devoted to brief summary. 



II. SERIES EXPANSION METHODS 

We begin . by b riefly summarizing the series expan- 
sion methodH3a We employ here the cluster expansion 
around a given strong-coupling spin singlet state. Let 
us explaiftrih e ii daa .taking the dimer expansion as an 
example.EjO e2lEa£a First the 2D Hamiltonian is divided 
into two parts: H = + Hi. The unperturbed Hamil- 
tonian Hq is composed of an assembly of isolated singlet 
dimers which are formed by the strong antiferromagnetic 
bonds. Namely, our starting configuration for the per- 
turbation has the disordered ground state with the spin 
gap. We then introduce the interaction term Hi among 
the independent dimers and observe how the physical 
quantities are changed by exploiting the power series ex- 
pansion with respect to Hi . The advantage of this cluster 
expansion is that we can combine analytical and numeri- 
cal techniques in a complementary way. For example, the 
computer can be utilized to systematicailM-a;enerate the 
higher order terms from the lower onesHJ'ESl We exploit 
the cluster expansions which may be most appropriate 
for each system with different structures. 

To discuss how the introduction of Hi perturbs the dis- 
ordered state with the spin gap and enhances the antifer- 
romagnetic correlation, we calculate the dispersion rela- 
tion i5(k) for 2D spin systems with various structures .1121 
This quantity is expanded as a power series in A as 
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aimn cos{lkx + mky)X^ , 



(1) 



where the wave number is denoted by k = [k^, ky) and 
the Brillouin zone for each model will be defined in the 
following section. We shall calculate the dispersion rela- 
tion up to the eighth order in A for the ladder-structure 
system and the fifth order for both the plaquette sys- 
tem and the mixed-spin system. To estimate the min- 
imum value of the dispersion relation in the first Bril- 
louin zone, we can also expand the spin gap A up to 
the same order. Since we are not able to analyze the 
critical phenomena only with the obtained power-series, 
further asymptotic analyses are necessary to discuss the 
phase transitions. To this end, we make use, of the Fade 
approximants and the differential methodsEJ to estimate 
the critical point for the phase transition, the disper- 
sion relation, etc. Especially, the critical point between 
the magnetically ordered and disordered phase is esti- 
mated not only by the ordinary Dlog Fade approximants 
but also by the biased Fade approximants. In the bi- 
ased method we assume that the phase transition in our 
2D quantum spin model should belong to the^niversal- 
ity class of the 3D classical Heisenberg modelE3 Namely, 
the critical value Ac for the perturbation parameter is de- 
termined by the formula A (Ac — A)'' with tiie known 
exponent — 0.71 around the transition point.c3 We also 
apply the first-order inhomogeneous differential method 
to the power-series to obtain the dispersion relation. It 
should be noted here that since higher-order coefficients 
in the series expansions are necessary to deduce the dis- 
persion relation in this method correctly, we might be 
sometimes left with wrong values at a certain wave num- 
ber after applying the asymptotic analysis. It is known 
that this type of pathology occasionally happens in these 
asymptotic approximations .EJ If we carefully discard this 
spurious behavior to find the correct one, these analyses 
provide a fairly good approximation in many cases, which 
will be explicitly shown in each case treated below. 



III. EXCITATION SPECTRUM AND PHASE 
DIAGRAM 

Let us now introduce the 2D antiferromagnetic quan- 
tum spin system defined by the Heisenberg Hamiltonian, 



H — Ho + Hi, 
Hq ^ Jl Si ■ Sj 

Hi — J2 ^ Si ■ Sj + J3 ^ Si ■ Sj 

(*,i)6G2 (»,i)GG3 



(2) 
(3) 

(4) 



where Ji, J2 and J3 denote the antiferromagnetic cou- 
pling constants, and Sj is the spin operator at the j-th 
site. To treat the mixed-spin systems as well as the lad- 
der and plaquette systems, the spin Sj is allowed to take 



different values at each cite. We denote the bonds 
for the non-perturbed Hamiltonian as Gi, while those 
for the perturbed parts as G2 and G3. By appropriately 
choosing the set of (Gi, G2, G3), we can deal with the 2D 
systems with various structures by the series expansion 
techniques. We treat below the case of A(= J2/J1) < 1 
and aA(= J3/J1) < 1(0 < a < 1). In the following, the 
excitation spectrum is analyzed to discuss the quantum 
phase transitions for the 2D antiferromagnetic spin sys- 
tems with ladder, plaquette and mixed-spin structures. 
We carry out the dimer expansion, the plaquette expan- 
sion and the mixed-spin cluster expansion. Starting with 
the above strong-coupling spin singlet states, we can per- 
form the cluster expansion with respect to A and aX. 



A. Dimer expansion for ladder-structure systems 

We first discuss a 2D spin system with the ladder struc- 
ture, which is shown schematically in Fig. ||, where the 
bold, the thin and the dashed lines represent the cou- 
pling constants 1, Al and olAl, respectively. It is noted 
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FIG. 1. The 2D s = 1/2 spin system with ladder structure. 
See the text as for the meaning of the bold, the thin, and the 
dashed lines. 

here that the spin ladder system (aL = lO^ was already 
studied in detail by the cluster expansion.Ej By changing 
Q!L, we can see how the isolated 2-leg ladder (q;l = 0) 
is changed to the 2D system. We calculate the energy 
for spin-triplet excitations by means of the dimer expan- 
sion up to the eighth order in Al for various values of 
Q!L. Note that the Brillouin zone is reduced to half of 
the original one because the dimer singlet is composed 
of two spins in the y-direction. By apolying the first- 
order inhomogeneous differential methodcj to the power 
series computed above, we obtain the spin-triplet excita- 
tion spectrum shown in Fig. ^ in the case of Al — 0.5. 
When q;l = 0, the system is reduced to the isolated 2- 
leg ladders with the inter-leg (intra-leg) coupling con- 
stant 1(Al = 0.5), which is knowja.to have the disordered 
ground state with the spin gap.o This gives rise to the 
flat dispersion between (tt, 0) and (7r,7r/2). The com- 
puted coefficients for the spin gap A = i?(7r, 0), in the 
series of Al are tabulated for some particular values of 
Q!L in Table ||. Note that for the isolated ladder case 
(ttL = 0), our results copespond to those obtained pre- 
viously by Oitmaa et alx^. The obtained spin gap with 
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FIG. 2. Plots of the spin-triplet excitation spectrum -E(k) 
along high-symmetry cuts through the Brillouin zone for the 
system with the couplings ql = 0.0,0.5,1.0 (shown in the 
figure from the top to the bottom at (tt, 0), respectively) when 
Al = 0.5. 

TABLE I. Series coefficients for the dimer expansion of the 
spin gap A = Eiji, 0) for the coupled-ladder system. 



n 


QL = 0.0 


QL = 0.2 


QL = 0.5 


QL = 1.0 





1.0000000 


1.0000000 


1.0000000 


1.0000000 


1 


-1.0000000 


-1.1000000 


-1.2500000 


-1.5000000 


2 


0.50000000 


0.38500000 


0.15625000 


-0.37500000 


3 


0.25000000 


0.16025000 


0.066406250 


0.031250000 


4 


-0.12500000 


-0.17002083 


-0.22737630 


-0.34635417 


5 


-0.27343750 


-0.22812954 


-0.26126692 


-0.88050673 


6 


-0.15332031 


-0.041639568 


0.018219038 


-0.16209751 


7 


0.24560547 


0.30931367 


0.36472170 


-0.20085681 


8 


0.48133850 


0.37176730 


0.34136063 


-0.63949210 



a fixed Al is shown in Fig. ||. It is seen that the spin 
gap decreases with the increase of the inter-ladder cou- 
pling Q!L and finally vanishes at which the quantum phase 
transition to the antiferromagnetically ordered state oc- 
curs. We wish to mention that the order in our cluster 
expansion is not high enough to deduce the accurate dis- 
persion for ttL close to the transition point within the 
first-order inhomogeneous differential approximation, as 
seen in Fig. I It thus seems difficult to deduce the crit- 
ical point Q!c correctly. However, as far as the critical 
value is concerned, we can use alternative analysis based 
on the Fade approximants, which provides a rather accu- 
rate estimate for cKc, by assuming A ^ [a.^ — aY near the 
critical point. By employing the latter analysis comple- 
mentarily around the critical point etc, we have obtained 
the corrected spin gap as a function of the inter-ladder 
coupling ckL, which is shown as the solid line in Fig. ^. 

We also show the phase diagram for the coupled- 
ladder system in Fig. |[ The solid line represents the 
phase boundary obtained by the biased [4/3] Fade ap- 
proximants for the spin gap, and the dashed line is 
the boundarji—determined previously by the staggered 
susceptibilityO. It is remarkable that the present resjaltria 
in fairly good agreement with the QMC simulations ,L2rL3 
and considerably improves the previous one especially in 
the region with small ckl. 

The above analysis may not be sufficient to discuss 
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FIG. 3. Spin gap A = i?(7r, 0) for the coupled-ladder sys- 
tem in Fig. |l|. The dots with error bars are the results 
obtained by the first-order differential method, whereas the 
solid lines denote the corrected values by applying the Fade 
approximants to the spin gap around the critical point com- 
plementarily. 
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FIG. 4. Phase diagram for the coupled-ladder system in 
Fig. ^. The solid line indicates the phase boundary ob- 
tained by biased [4/3] Fade approximants for the spin gap 
while rthe dashed line have been obtained in our previous 
paper Ji3 The-i^id circles represent the results of the QMC 
simulations.! 



the 2D ladder-structure systems generically, because the 
parameter region treated by the cluster expansion is re- 
stricted. To make a complementary analysis, we next 
regard the present system as the coupled-dimer chains, 
for which the bold, the thin and the dashed lines in Fig. 
|l| represent the coupling constants 1, cidAd and Ad- We 
carry out the similar calculation up to the eighth order 
in aD, and list the resulting power series for several val- 
ues ttD in Table ||. Applying the Fade approximants 
to the computed spin gap, we obtain the phase diagram 
shown in Fig. ^. In this figure, the solid line represents 
the phase boundary obtained by the biased [4/3] Fade 
approximants, and the dashed line is the-pne obtained 
previously by the staggered susceptibility.Ea We find that 
these two boundaries are in fairly good agreement with 
each other, and fuilhermore consistent with those of the 
QMC simulationsEj'O (the solid circles in Fig. ^). By 
these comparisons, we can say that our cluster expan- 
sion approach gives quite accurate results for the phase 
diagram. 
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TABLE II. Series coefficients for the dimer expansion of 
the spin gap A — E{tt, 0) for the coupled-dimer-chain system. 



n 


ao = 0.0 


od = 0.2 


an = 0.5 


QD = 1.0 





1.0000000 


1.0000000 


1.0000000 


1.0000000 


1 


-0.50000000 


-0.70000000 


-1.0000000 


-1.5000000 


2 


-0.37500000 


-0.45500000 


-0.50000000 


-0.37500000 


3 


0.031250000 


0.063250000 


0.062500000 


0.031250000 


4 


-0.013020833 


-0.030220833 


-0.067708333 


-0.34635417 


5 


-0.061930339 


-0.10737895 


-0.25710720 


-0.88050673 


6 


0.010735971 


0.040159654 


0.039823179 


-0.16209751 


7 


0.0030713964 


-0.017981720 


-0.12801535 


-0.20085681 


8 


-0.031547664 


-0.077754777 


-0.28991249 


-0.63949210 




CCD 

FIG. 5. Phase diagram for the coupled-dimer-chain sys- 
tem in Fig. |l| The solid line indicates the phase boundary 
obtained by the biased [4/3] Fade approximants. We also 
show the phase boundary which was determined by applying 
Dlog [4/3] Fade approximapts to the fourth-order series for 
the staggered susceptibiUtafJij The solid circles represent the 
QMC simulation results.L3'L3 



B. Plaquette expansion 

In the following, we consider the plaquette-structure 
systems. Introducing the spin systems with two kinds 
of the plaquette structures, we discuss the quantum 
phase transitions between the ordered and the disordered 
states. We note here that series expansion studies on 
plaquette systgrns, lifive been done extensively by several 
groups so far,Q'Era'Ej which we shall also compare with 
our results in some special cases. 



1. plaquettes on a square lattice 

First, we treat the plaquettes on a square lattice shown 
in Fig. |[ The starting Hamiltonian Hq is composed of 
the isolated plaquettes, whose ground state is spin sin- 
glet with the excitation gap A = 1. We study how the 
antiferromagnetic correlation develops in the presence of 
the inter-plaquette interaction A and aX. In the previous 
paper jl3 we calculated the staggered susceptibility up to 
the fourth order and determined the phase boundary be- 
tween the magnetically ordered and the disordered states 
(see the dashed line in Fig. ph. We here calculate the dis- 
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FIG. 6. Flaquette structure on a square lattice. The bold, 
thin, and dashed lines denote the coupling constants 1, A, and 
a A among s = 1/2 spins 

TABLE III. Series coefficients for the plaquette expansion 
of the spin gap A = E{0, 0) for the plaquette system on a 
square lattice. 



n 


a = 0.0 


a = 0.2 


a = 0.5 


a = 1.0 





1.0000000 


1.0000000 


1.0000000 


1.0000000 


1 


-0.66666667 


-0.80000000 


-1.0000000 


-1.3333333 


2 


0.019675926 


-0.068425926 


-0.19762731 


-0.40509259 


3 


0.064935378 


0.016056713 


-0.081316913 


-0.28178048 


4 


0.043061549 


0.019385674 


-0.038615359 


-0.20391542 


5 


0.039539538 


0.030819425 


-0.019630017 


-0.23535878 



persion for spin excitations up to the fifth order, and list 
the obtained series of the spin gap A = E{0, 0) for sev- 
eral values of a in Table |l|. Using the first-order inho- 
mogeneous differential methods, we obtain the dispersion 
relation shown in Fig. ^. Note that the Brillouin zone is 



1 .75 
1 . 5 
1 .25 

1 

.75 
0.5 
0.25 



P-1 



(0,0) (Jt/2,0) (7U/2,7U/2) (0,0) 

k 

FIG. 7. Flots of the spin-triplet dispersion ECk) along 
high-symmetry cuts through the Brillouin zone for the pla- 
quette structure system with the couplings a = 0.0,0.5,1.0 
(shown in the figure from the top to the bottom at (0,0), 
respectively) in the case A — 0.5. 

reduced to a quarter of the original one due to the plaque- 
tte structure. We also show the spin gap A — E{0, 0) as 
a function of a in Fig. ||. As the inter-plaquette coupling 
a is increased, the antiferromagnetic correlation grows 
up, which causes the decrease in the excitation gap, and 
finally induces the quantum phase transition to the an- 
tiferromagnetic with the vanishing spin gap. In the case 
{a, A) = (0, 1), our model is reduced to the independent 
isotropic two-leg ladders for which the spin gap A = 0.43 
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FIG. 8. Spin gap A — E{0, 0) for the plaquette system in 
Fig. ^. The dots denote the results computed by the first 
order differential method, whereas the solid line represents 
those corrected with the Fade analysis around the critical 
point. 



is obtained. This value is slightly small compared—with 
A = 0.504 (density matrixr*enormalization group)u and 
0.5028 (dimer expansion) ,£21 which implies that higher- 
order cluster expansions may be necessary to obtain more 
accurate values of the spin gap for the plaquette system. 
In contrast, it is shown below that the phase diagram 
can be obtained with much higher accuracy. By applying 
Pade approximants to the power series of the spin gap, 
we obtain the phase diagram in Fig. ^ It is remarkable 




a 

FIG. 9. Phase diagram for the plaquette system on a square 
lattice. The solid Hne indicates the phase boundary obtained 
by the Dlog [2/1] Fade approximants and the dashed lincris 
the one obtained previously by the staggered susceptibilityJij 
The solid circles denote the QMC simulation results.EiJ 

that the phase boundary given in this paper considerably 
improves the previous one in the small a regime, which 
can be confirmed by the result-ef the QMC simulations 
(the dot shown in the figure). Ell We note here that for 
the special case of a = 1, sirpalar results were praviously 
reported by Fukumoto et alu and Wcihong et alB 



2. plaquettes on a 1/5 depleted square lattice 

We next deal with the plaquette system shown in Fig. 
p^ , which may be regarded as a 1/5 depleted square 
lattice,ErEl by extending the work done by Gelfand et al.tl 
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FIG. 10. 2D spin system on a 1/5 depleted square lattice. 
The bold, the thin, and the dashed lines indicate the coupling 
constants 1, A, and qA among s = 1/2 spins. For simplicity, 
we define the lattice constant as the distance between plaque- 
ttes. 

TABLE IV. Series coefficients for the plaquette expansion 
of the spin gap A = E{tt, tt) for the plaquette system on a 
1/5 depleted square lattice. 



n 


a = 0.0 


a = 0.2 


Q = 0.5 


a = 1.0 





1.0000000 


1.0000000 


1.0000000 


1.0000000 


1 


-0.33333333 


-0.40000000 


-0.50000000 


-0.66666667 


2 


-0.10590278 


-0.13236111 


-0.18793403 


-0.32291667 


3 


0.032959989 


0.024330922 


0.0092901536 


-0.0081862461 


4 


0.024698692 


0.022081937 


0.020367093 


0.033841824 


5 


-0.0070637727 


-0.0088223796 


-0.014159237 


-0.044030516 



This system is-^bg-cjeiisidered to be made out of the pla- 
quette chains Jj'Eaoo since the model with a = and 
finite A is reduced to the isolated plaquette chains. The 
results for the cluster expansion of the spin gap up to the 
fifth order are tabulated in Table IV for several values of 



a. The resulting value of A deduced by the Pade anal- 
ysis is shown in Fig. ^ as a function of a. In the case 
of a = 0, the model is reduced to the isolated plaquette 
chains with the spin gap. In this case, by applying the 
differential methods to the power series, the spin gap is 
estimated as A = 0.607 ± 0.001 for A = 1, which is in 
good agreement with— the result of the exact diagonal- 
ization A = 0.6086JjO To observe the phase transition 
when the couplings a between the plaquette chains in- 
creased, the phase diagram is shown in Fig. Here, 
the phase boundary (solid line) is determined by apply- 
ing the biased [2/3] Pade approximants to the spin gap. 
We find that this line is quite consistent with the previ- 
ous resultstj shown as the dashed line. The fact that the 
two lines evaluated for different quantities in different or- 
ders produce a quite similar behavior confirms that the 
obtained boundary is indeed reliable although our cal- 
culation is restricted to the lower-order expansions. We 
also note that the results already obtained by QMCq and 
also by the plaquette expansionQ in the case of a = 1 are 



5 
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FIG. 11. Spin gap A = E{-k,-k) for the plaquette system 
in Fig. ^. The dots are the data computed by the first order 
differential method, and the solid line denotes the corrected 
one by combining the Fade analysis around the region with 
the small spin gap. 
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FIG. 12. Phase diagram for the plaquette system in Fig. 
p^ . The solid line denotes the phase boundary obtained by 
the biased [2/3] Fade approximants while the dashed line is 
the previous resultsc-l obtained by the staggered susceptibilitsi 
The solid circles denote the results of the QMC simulation.u 



in good agreement with the present one. 



C. Mixed-spin cluster expansion 

Let us now turn to another interesting 2D system com- 
posed of two kind of different spins, which has attracted 
much attention recently. In this mixed-spin system, the 
topological nature of spins is important for the system 
to generate the spin gap or induce an antiferromagnetic 
long-range order. jIh this subsection, we extend the pre- 
vious calculationslla to those of the excited states, and 
quantitatively discuss the phase transition in 2D mixed- 
spin systems. We will clarify that the arrangement of dif- 
ferent spins affects the nature of the quantum phase tran- 
sitions from the spin-gap phase to the antiferromagnetic 
phase. We shall also check that our series expansion ap- 
proach correctly captures the spin structure though our 
calculation is based on the lower-order perturbations. 

We deal with two typical systems composed of s = 1 /2 
and 1, as displayed in Figs. |l^ and |lJ.Ej 



FIG. 13. 2D spin system with the columnar-type 
mixed-spin structure with s = 1,1/2. We denote the cou- 
phng constants 1, A, and aX. by the bold, the thin, and the 
dashed lines. 
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FIG. 14. 2D spin system with the 
structure. The meanings of the bold, the thin, and the dashed 
lines are the same as those in Fig. hsl 



1. columnar-type mixed-spin system 

We begin with the columnar-type mixed-spin system, 
for which the mixed-spin chains are stacked uniformly in 
a vertical direction (Fig. |l^). Starting from the mixed- 
spin clusters of 1/2 o 1 o 1/2, we perform the series ex- 
pansion with respect to A. Note that the Brillouin zone 
is reduced to a third of the original one since the mixed- 
spin cluster is composed of three spins in the x-direction. 
We list the power series obtained up to the fifth order for 
the excitation spectrum in Table It is noted that in 
the case of the isolated mixed-spin chain (a = 0), these 
coefficients are the same as those for the plaquette chain, 
(see Fig. |l^) and thus the isotropic mixed-spin chain 
with A = 1 has the same spin gap A = 0.607. We can 
indeed prove that the mixed-spin chain is identical to the 
plaquette chain as far as the ground state and the low- 
energy elementary excitation are concerned. Using the 
first-order inhomogeneous differential methods, we ob- 
tain the dispersion relation shown in Fig. ^ We recall 

TABLE V. Series coefficients for the mixed-spin cluster ex- 
pansion of the spin gap A = _E(7r/3,7r) for the 2D columnar 
mixed-spin system. 



a = 0.0 



a = 0.2 



a = 0.5 



a = 1.0 



1.0000000 
-0.33333333 
-0.10590278 
0.032959989 
0.024698692 
-0.0070637727 



1.0000000 
-0.73333333 
-0.20868056 
-0.014472881 
0.013380005 
-0.046571788 



1.0000000 
-1.3333333 
-0.24826389 
-0.090593252 
-0.14493772 
-0.44465163 



1.0000000 
-2.3333333 
-0.008680556 
0.50585737 
-1.8264592 
-5.2417951 
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(0,0) 



FIG. 15. Plots ol the spin-triplet dispersion relation 
ii'(k) along high-symmetry cuts through the Brillouin zone 
for the columnar mixed-spin system with the couplings 
a = 0.0, 0.3, 0.6 (shown in the figure from the top to the 
bottom at (7r/3,7r), respectively) in the case A = 0.5. 



that the mixed-spin system in the case of a = is re- 
duced to the mixed-spin chain with the spin gap defined 
at the wave number k = (tt/S, tt). Increasing the inter- 
chain coupling a, we can see that the spin gap decreases 
as the magnetic correlation grows up, and finally van- 
ishes at which the phase transition to the magnetically 
ordered phase takes place. In Fig. |l^, the phase bound- 



1.2 




a 

FIG. 16. The left (right) solid line represents the phase 
boundary for the mixed-spin system in Fig. |l^ ( Fig. 
The results obtained from the ground-state susceptibilityEj 
are shown as the dashed lines. 

ary is shown by the left solid line, which is obtained by 
the biased [2/3] Pade approximants for the spin-gap. In 
comparison, we also display the previous resultEj by the 
left dashed line, which was determined by the staggered 
susceptibility in the fourth-order expansion. As has been 
the case for the plaquette systems, we can see again that 
the phase boundaries which were determined via the dif- 
ferent physical quantities are consistent with each other. 
This demonstrates that the reliable phase boundary is 
established by the present analysis. 



2. diagonal-type mixed-spin system 

We next discuss the diagonal-type mixed-spin system 
shown in Fig. |l^, for which the mixed-spin chains are 
stacked diagonally. According to this structure, the 
shape of the Brillouin zone for the diagonal system is 
quite different from those for the columnar one as shown 
in Fig. |lj. The definition of the coupling constants is 




FIG. 17. First Brillouin zone for the 2D system with the 
mixed-spin structure. 



the same as that in Fig. 13, The mixed-spin cluster ex- 



pansion for the excited states up to the fifth order with 
the asymptotic analysis yields the phase diagram and 
the dispersion relations shown in Figs. |l^ and |l^, re- 
spectively. In Fig. |l^, the right solid line represents the 
phase boundary determined by the biased [3/2] Fade ap- 
proximants for the spin gap . The resulting series for 



2 

1.75 
1.5 

1.25 



0.75 



0.5 
0.25 




FIG. 18. Plot of the spin-triplet dispersion relation E(k.) 
for the diagonal mixed-spin system with the coupling param- 
eters a — 0.0, 0.3, 0.6 (shown in figure from the top to the 
bottom at (7r/3,7r/3), respectively) when A — 0.5. 



some particular values of a are tabulated in Table VT 
For a = 0, the system correctly reproduces an assembly 
of independent mixed-spin chains which have the disor- 
dered ground state. Increasing a, the correlation among 
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TABLE VI. Series coefficients for the mixed-spin cluster 
expansion of the spin gap A = i?(7r/3, 7r/3) for the 2D diago- 
nal mixed-spin system. 



: 0.0 



0.2 



: 0.5 



1.0 



1.0000000 
-0.33333333 
-0.10590278 
0.032959989 
0.024698692 
-0.0070637727 



1.0000000 
-0.66666667 
-0.22430556 
0.025971997 
0.0053494949 
-0.0082040994 



1.0000000 
-1.1666667 
-0.42925347 
-0.066280914 
-0.11649536 
0.020462458 



1.0000000 
-2.0000000 
-0.84375000 
-0.41995001 
-0.96473863 
-0.0043513404 



the mixed-spin chains grows up, and the quantum phase 
transition occurs. Especially, in the case of the mixed- 
spin chains with the isotropic bonds (A = 1), the phase 
transition to the ordered state occurs at the critical value 
ttc = 0.21. We note that as well as the case of the colum- 
nar case, this line is consistent with the phase hciundary 
determined from the ground-state susceptibilityjlj shown 
as the dashed line in Fig. which may ensure that 
phase diagrams for both of the two distinct mixed-spin 
systems are determined in rather high accuracy. 



IV. SUMMARY 

We have performed the systematic cluster expansion 
to study the two dimensional quantum spin systems with 
modulated lattice as well as spin structures. By apply- 
ing the asymptotic analysis to the obtained series, we 
have calculated the excitation spectrum and have dis- 
cussed the quantum phase transitions. We have thus con- 
structed the phase diagram which improves the previous 
one obtained via the staggered susceptibility. In particu- 
lar, we find that the present results for the systems with 
the ladder and plaquette structures are in fairly good 
agreement with the results of the QMC simulations. We 
have further studied the critical phenomena for the spin 
systems with modulated S* = 1,1/2 structure. By care- 
ful study on two types of slightly different systems with 
mixed-spins, it has been clarified that the topology of the 
spin arrangement plays an important role to stabilize the 
spin-gap phase. 
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